Integrable model for interacting electrons in metallic grains 
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We find an integrable generalization of the BCS model with non-uniform Coulomb and pairing 
interaction. The Hamiltonian is integrable by construction since it is a functional of commut- 
ing operators; these operators, which therefore are constants of motion of the model, contain the 
anisotropic Gaudin Hamiltonians. The exact solution is obtained diagonalizing them by means of 
Bethe Ansatz. Uniform pairing and Coulomb interaction are obtained as the "isotropic limit" of 
the Gaudin Hamiltonians. We discuss possible applications of this model to a single grain and to a 
system of few interacting grains. 
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Introduction and summary of the results. Progress in 
nanotechnology has opened up theoretical investigations 
on the behavior of disordered interacting systems of small 
size Q. Recently, the I-V characteristic measurements 
of Ralph, Black and Tinkham Q on small Al dots stim- 
ulated the theoretical debate on how to characterize the 
physical properties of small metallic grains, such as su- 
perconductivity and ferromagnetism (HQl. Due to the 
chaoticity of the single particle wave functions , the 
Hamiltonian of these systems reads 



J E4^-'^v + 0(6E 2 /E T ). (1) 



(Here and in the following, sums over spins a, a' are im- 
plied). The quantum numbers i, a label a shell of doubly 
degenerate single particle energy levels with energy ej and 
annihilation operator Ci CT ; :— c\ a Ci a \ S a , a = x,y,z, 
are 2x2 spin matrices; SE is the average level spacing, 
and Et the Thouless energy. The universal part of the 
Hamiltonian ([l]) (namely the first four terms) describes 
the pairing attraction, the electrostatic interaction and 
the ferromagnetic instability, respectively. 
The superconducting fluctuations [0 can be taken into 
account by employing the BCS model |3]|| (namely tak- 
ing the first two terms in Eq. (|l|)). Richardson and Sher- 
man (RS) 0| constructed the exact solution of the BCS 
model by a procedure close in spirit to the coordinate 
Bethe Ansatz (BA). The knowledge of the exact eigen- 
states and eigenvalues of the BCS model has been crucial 
to establish physically relevant observables [|| . The inte- 
grability of the model has been proved |) 10 1 to be deeply 
related to the integrability of the isotropic Gaudin mag- 
net |yj: the BCS model can be expressed as a certain 
combination (see Eq. @ below) of its integrals of motion, 
which contain Gaudin Hamiltonians. Relations with con- 
formal field theory and disordered vertex models were in- 
vestigated in Refs. |l2]Jl3[ ]. 

Many properties of metallic grains in a normal state (neg- 



ligible superconducting fluctuations) can be described by 
the orthodox model jj],[l4| (i.e. taking the first and the 
third term of the Hamiltonian (Q)). This arises by as- 
suming uniform Coulomb interaction. 
Magnetic phenomena like the mesoscopic Stoner instabil- 
ity |H] can be studied by means of the exchange contri- 
bution to the Hamiltonian (the fourth term in Eq. (|I|)). 
The terms proportional to 5E 2 /Et correspond to non 
uniform Coulomb interaction Although they lose 

importance with the increasing conductance of the sys- 
tem, these corrections gain physical relevance due to the 
typically low relaxation rate of the excitations in a small 
dot. In fact the corrections to the orthodox model induce 
"fluctuations" which can explain how non-equilibrium ex- 
citations decay in the dot |j^Jr7|j . This results in the 
formation of clusters of resonance peaks in the tunneling 
spectroscopy experiments 

In this paper we present an integrable generalization of 
the BCS Hamiltonian with non-uniform pairing coupling 
gij and solve it exactly. Besides the non-uniform pairing, 
the Hamiltonian contains a non-uniform Coulomb inter- 
action Uij ; and fTy are fixed according to Eqs. (||) . Wc 
shall see that the inclusion of certain 0(SE 2 / Et) terms 
leads to our integrable model. The integrable Hamilto- 
nian we solve is 



crcr'^jcr' I 3 



(2) 



where the couplings are 

gij = -qK(si - e J )/sinhg(u l - uj) , 



A + qK(ei — £j) cothq(ui — Uj), 



9jj = -P. 



3 ! 



Wjj = A + (3j, 



(3) 



where 20 j = -qKj2 i ^ j (s i -e 3 ) cothg(w l -Uj) + C. Q 
For generic choices of ]3j, the single particle energies Ej 
must be shifted by 2(5 j + 4 Ylijtj Uij in order to have in- 
tegrability. The parameters A, K, and C are arbitrary 



1 



real constants, while q can be real or imaginary. The BCS 
Hamiltonian, including a tunable capacitive coupling can 
be obtained from (|J) in the isotropic limit q — > 0. Non- 
uniform coupling constants are obtained for generic q, 
and Uj being monotonic functions of Sj. For real q, the 
arising g^j can be made nearly uniform for levels within 
an energy cutoff Ed, and exponentially suppressed oth- 
erwise; correspondingly, [Zy can be made nearly uniform 
(as specified in Eqs.(p0|) below). 

The proof of the integrability of the Hamiltonian (||) pro- 
ceeds along the two following steps, i) First we note 
the factorization of the eigenstates of the Hamil- 
tonian (H): I*) = \^n) <8> |$m) with eigenvalue E = 
E^ + Em', where \^>n) is the eigenstate of the Hamilto- 
nian Hpf projected on the subspace with N time-reversed 
pairs; |$m) is the Fock state projected on the blocked 
M singly occupied levels. The solution of the corre- 
sponding Hamiltonian Hm is easily obtained p(J as: 
H M |$m) = Ei Ei + E« Ua - JS(S+ 1)] |$ M >. ii? Then 
Hamiltonian (0) is integrable if and only if iijv is inte- 
grable. The Hamiltonian Hn is obtained by inverting 
the procedure presented in Ref. [ju^: First, we modify 
the constants of motion (of the BCS model) to commut- 
ing operators containing the anisotropic Gaudin models 
(the isotropic ones being considered in Jl(]]); then we de- 
fine the Hamiltonian in terms of these operators (Hn is 
therefore integrable by construction). We discuss some 
choices of {uj}, K, and A leading to physically relevant 
Hamiltonians. The exact solution of Hn is found by di- 
agonalizing the integrals of motion through BA |ll[] . The 
exact eigenf unctions and eigenvalues ^n, En are 



a—l j—1 J 

n n 

En = qK Sj coth q((J a — u 

j=l a=l 



AN 2 



(4) 
(5) 



solution is presented. This will complete the study of 
the integrability of the Hamiltonian (g). Finally, we will 
explain how our model can be applied to describe single 
as well many interacting grains. 

Integrability. The BCS Hamiltonian can be written in 
terms of the spin- 1/2 realization of su(2): 

Hbcs = 2£ i s 3 -aY, s t s k' where 

j j,k 

Sj ■= c n c jh S+ := (Sr)t = ct^, 
SI~\{c)^ + c] iCji -l), (7) 

obeying [S],S±]= ±6 jk S£, [S+,S^] = 25 jk S%. Its 
constants of motion are written in terms of isotropic 
Gaudin Hamiltonians S, 



k = l 3 
k^3 



(8) 



The Tj mutually commute and we have [fj , ffc] = 
[H , fj] — for all i,j e {1, - - - , 0}, because the BCS 
Hamiltonian can be written in terms of the fj only: 



H 



BCS 



Z_^,^3'3 

3 



(9) 



3,k 



Our approach is now to modify the integrals of motion 
(||) and then to construct an integrable BCS-like model 
(which turns out to be characterized by a non-uniform 
pairing) following formula (Q): 



H 



N 



E 2 * 



.4 



\ " 



3,k 



TjTk + const. 



(10) 



The ansatz for the modified integrals Tj is 



|0) is the electronic vacuum state and fl is the number of 
levels. The quantities uj a fulfill the equations 



Sj + a j 



o _ \ ^ a na no 



Z— I J* 3 

k = l 



(11) 



— - ^Tq coth q(cj a - 



u,) 



1=1 



N 



1 cothg(w a - wp) = ,a = 1,..., N (6) 



0=1 



Our results can be applied to describe a system of TV 
grains, since their Hamiltonian can be written (after a 
suitable relabeling of the levels) in the form (g). For 
distinct grains gjk describe Cooper pair tunneling, and 
Ujk the inter-grain Coulomb interaction. We require gij 
to decay both with inter-grain distance and level separa- 
tion. This can be fulfilled with Uj fixed by Eqs.(22), (23). 

The present paper is laid out as follows. First we dis- 
cuss the integrability of the model Hn', then its exact 



where the operators Sj are anisotropic Gaudin Hamilto- 
nians (the isotropic case corresponding to w 
wfj). These operators mutually commute if 



x _ ...v _ 



W ij W jk + W ji W ik = W ikWj k , 



"jll 



(12) 
(13) 



where ( |l2| ) emerges from imposing [Si,Sj] = pl| . 
The other condition arises from [Sf,Sj] + [S», £?] = 0. 
We furthermore postulate particle number conservation, 

n 

which in the spin picture means [^2 Sf,^] = for all 

i=l 

j € {1, . . . , n}, leading to another condition 



Wij = w tJ 



y Eq^l3) 



-II),- 



-Wj 



(14) 



2 



The last equation reduces the anisotropy to the XXZ- 
type and Eqs. (12) finally become 



WijVjk + WjiVik = WikWjk , Vij := i 
The solution of Eqs. @ (see Ref. Q) is 
Vjk = qK coth q(uj - Ufe), 



(15) 



sinhg(uj 



life) 



(16) 



where tij are arbitrary complex parameters such that 
Vjk,Wjk are real. The transition from hyperbolic to 
trigonometric functions in the solution Jl6|), is gained 
through the choice q = i, with real K , uj. The cubic and 
quartic terms in S? (obtained from formula (|l0|)) vanish 
for the antisymmetry of Vjk- We finally obtain: 



j 



J2 mSpk + 4 £ U jk S*S z k , (17) 

j,k j,k 



where the couplings are given in Eqs. (Q). Up to a con- 
stant, the Hamiltonian (0) (projected on doubly occu- 
pied states) is recovered writing back the spin operators 
in terms of creation and annihilation operators. 

Exact solution. The exact solution of the anisotropic 
Gaudin model for my and My fixed by Eqs. ( |l6| ) was ob- 
tained in Ref. pT| . The same procedure can be applied 
to diagonalizc Tj. The cigcnfunctions of Tj defined in 
Eq.(pl|) are written in the form 



l*i> 



c 0i>-- 

3i<-"<3M 

£' e ^ 

3i<-..<jm-i 



,j M )St...St t \0) (18) 
..,jM-i)S+...S+ u _ x S}\0). 



The vacuum |0) corresponds to |j, . . . , J.) ; the prime on 
the sums means the indices run in the range {1, . . . , Q}/ 
{j}. Imposing that \tyj) is an eigenstate of Tj we find 
a set of equations which c({jj}) and e({jj}) must ful- 
fill. For a suitable change of variables we find that these 
conditions are transformed in Eqs.(Q). The quantities Tj 
have the following eigenvalues: 



1 1 n ' 

Tc h * = 9E q coth ?( u a 



Z=l 



AT 

M;) — g coth g(w 

a=l 



(19) 

■3 w q) ■ 

t). The 



The parameters w a are determined by Eq.(|^ 
eigenvalues of Hm immediately follow from formula (10) 



Together with the eigenfunctions they are given in 
Eqs.(@),(|). 

Single grain. We discuss how our results can be ap- 
plied to describe the physics of a single grain. 



The isotropic limit q — > of Eqs. (^|) gives the BCS Hamil- 
tonian plus a tunable capacitive coupling A+g, with K = 



9 /Ed, Pi 



-9, 



-e 3 /(E D 0(|e, -E F \ -E D )), 



where Ep is the Fermi level, and is the Heaviside func- 
tion (0(x) = 1 if x < 0, 0(x) = if x > 0), setting sharp 
cutoffs at the Debye energy (2l]]; the diagonal elements 
Ujj and gjj can be independently set to arbitrary values 
(since they would renormalize Sj). Choosing A = —g 
gives the "pure" BCS model. In this limit, the eigen- 
states and eigenvalues Eqs.(0),(||) coincide with those of 
the BCS model and Eqs. (g) reduces to the RS equa- 
tions @. 

We now discuss the case corresponding to q = 1: 



K = g/ Ed , 0i = -9 , A > (g/E D ) max{fr 3 - e k } 

0,k 



-*i I Ed 



(20) 



We can identify three regimes depending on the value of 
Ed- i) Ed < SE, gij is nearly zero, while E/y ~ A — g; 
ii) Ed ~ SE, the pairing interaction decays on the scale 
Ed ~ SE, while f/y is slowly modulated by the energy 
separation; Hi) Ed > maxi.j(£i — Sj) both <?y and Uij 
are nearly uniform. 

Application to many interacting grains. We now dis- 
cuss applications of the model (||) to interacting dots. 
The Hamiltonian (g) can be reinterpreted as follows: the 
set I — {1, ...,fl} can be split into the (disjoint) sets 
I a , a = 1 . . . Af containing the levels of the a-th grain: 
I = U a 4; n = Y^=i&a, where fl a = \I a \. Thus the 
Hamiltonian Hm is equivalent to the following one: 



(21) 



AT 



EV- a,h f f 

a, 6—1 i a 
AT 

a, 6—1 z a ,i& 

where i a = 1, . . . fi a label the elements of I a and c a ,i a a 
annihilates an electron with spin a in the i a -th level of 
a-th grain. For a ^ b, g( a ' b ) describe tunneling of Cooper 
pairs; in terms of (|^), g^'^J = gij, where i is the z Q -th 
element of I a , and j the jh-th element of 7^; U^ "^ de- 
scribe a Coulomb-like coupling between grains a and 6, 
and is written in terms of Uij analogously to g( a ' b K Cou- 
plings g( a ' a ) and £/( a,a - ) describe pairing and Coulomb 
intra-grain interactions respectively. We fix the couplings 
as in ( |20|) with the exception that 



= 5> a - Sj/E D , when j e I a . 



Now we impose 



$ a > max {(ej - e k ) /E D } 



(22) 



(23) 
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to make the tunneling amplitude exponentially sup- 
pressed with the spatial distance between the grains. The 
pairing interaction is nearly uniform for levels within Ed 
in the same grain. The intra-grain Coulomb interaction 
is also nearly uniform Ujk — A, while the inter-grain 
Coulomb interaction is modulated by the corresponding 
energy separation. 

Conclusions. We found a class of integrable Hamilto- 
nians, which are a generalization of the BCS Hamiltonian 
characterized by non- uniform coupling constants. To our 
knowledge this is the first exact solution for non-uniform 
pairing interaction. The strategy we have adopted con- 
sists in generalizing the procedure of Ref. JlO[ , namely 
constructing the Hamiltonian of the system in terms of 
anisotropic Gaudin Hamiltonians. By means of the in- 
tegrability and the exact solvability of the latter we ob- 
tain the integrability and the exact solution of the model 
Eq. (|J). In this sense, our procedure is close in spirit 
to the quantum inverse scattering method |22[ . 
The isotropic limit q — > of the Gaudin Hamiltonians 
corresponds to uniform couplings. For arbitrary A, the 
Hamiltonian is the sum of the BCS and the orthodox 
model. For A — g the BCS Hamiltonian is obtained; the 
same isotropic limit of the exact solution Eqs. ([|)- @ 
coincides with the RS solution. 

This class of models might be useful for applications to 
the physics of metallic grains. The non- uniformity |25| ] 
of the coupling constant s (H ) corresponds to include cer- 
tain 0(SE 2 /Et) terms |l5[ in the Hamiltonian (0). In 
fact, we recover the "fluctuations" of the Coulomb inter- 
action of the Ref. jLj| identifying 8U h = — Uij< . The 
integrable model presented here might be applied as a 
starting point for suitable perturbation schemes leading 
to the explanation of the tunneling phenomena. 

The present model can be applied to systems of few 
interacting dots, since our capacitive-like inter-grain in- 
teraction does not decay with spatial distance. 

In a recent paper Ref. p4j| a non-uniform coupling for 
bosonic systems was studied. The Hamiltonian was con- 
structed from the bosonic analog of formulas jl0|), (11), 
where the S a are generators of su(l, 1) (instead of su(2)). 
This algebraic difference does not affect the equations 
which Wij, Vij have to fulfill to ensure the commutativ- 
ity of the (bosonic) tj. The coupling constants of this 
bosonic model can be obtained in the isotropic limit of 
our Eqs. (|J) with Uj oc ej and A = 0. This shows that 
the bosonic Hamiltonian in Ref. |^4j can be obtained by 
the limit q — > of anisotropic su(l, 1) Gaudin models. 
Work is in progress along this direction. 
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